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ABSTRACT 


An  experimental  study  of  waterfloods  was  conducted 
on  a  Berea  sandstone  core  to  evaluate  the  effects  of  initial 
gas  saturation  on  the  oil  recovery  and  to  study  the  mechanism 
of  the  waterflood  in  the  presence  of  the  gas  phase. 

Gas  saturations  ranging  from  zero  to  33  percent  were 
created  by  depletion  drive  performance  and  the  waterfloods 
were  conducted  at  five  pounds  pressure  differential  across 
the  core  32  centimeters  long  and  5  centimeters  in  diameter. 

Results  show  that  the  total  recovery  (depletion  drive 
plus  water  drive)  increases  steadily  from  42.75  percent  to 
64.65  percent  of  initial  oil  in  place  with  an  increase  of 
initial  gas  saturation  from  zero  to  33.4  percent. 

Also  it  is  found  that  the  presence  of  initial  gas 
saturation  up  to  approximately  10  percent  does  not  basically 
change  the  mechanism  of  displacement. 
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INTRODUCTION 


At  a  certain  time  in  the  life  of  a  reservoir, 
production  by  primary  methods  reaches  an  economic  limit 
and  it  becomes  necessary  to  rely  on  some  secondary  recovery 
techniques.  Among  these  water  flooding  has  been  given  a 
great  deal  of  attention.  Laboratory  and  field  experiments 
have  been  conducted  to  probe  into  the  problems  concerning 
water  flood.  All  possible  factors  affecting  the  performance, 
such  as  pressure  gradient,  velocity  of  water  flood,  initial 
fluid  saturations,  viscosities  of  the  fluids,  capillary 
forces  of  the  system,  wettability  of  the  rock  and  saturation 
history  have  been  investigated. 

A  number  of  investigators  have  tried  to  evaluate  the 
effects  caused  by  the  presence  of  free  gas  in  the  system 
before  water  flooding.  Many  of  the  experiments  were  con¬ 
ducted  on  sandstones  and  it  has  been  reported  that,  in  general, 
a  lower  residual  oil  saturation  was  obtained  by  water  flooding 
in  the  presence  of  free  gas  saturation  than  by  flooding  the 
same  system  in  the  absence  of  the  gas  phase.  The  degree  of 
improvement  in  recovery  varied  widely  depending  upon  the 
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system  and  the  conditions  of  the  test.  The  increased  oil 
recovery  obtained  because  of  the  presence  of  gas  during  the 
water  flood  has  been  attributed  to  changes  in  the  physical 
characteristics  of  oil,  selective  plugging  action  of  the 
gas,  inclusion  of  the  oil  mist  in  the  free  gas  phase, 
additional  sweeping  action  of  the  free  gas  and  simple 
replacement  of  the  oil  by  the  gas  phase.  It  is  visualized 
that  the  gas  will  exist  as  a  bubble  inside  the  discontinuous 
residual  oil  with  the  size  of  the  oil  bubble  substantially 
unchanged  by  the  presence  of  this  gas  bubble. (see  Figure  1) . 

It  is  also  thought  that  the  presence  of  a  free  gas 
saturation  inside  the  oil  will  reduce  the  relative  permeabil¬ 
ity  to  oil  which  will  exist  at  any  particular  water  satura¬ 
tion.  The  reduction  in  permeability  will  be  caused  by  two 
factors  -  the  additional  gas-oil  interface  and  the  reduction 
of  area  available  for  oil  flow  in  the  pores  containing  gas. 

Difficulties  are  encountered  in  attempting  to  arrive 
at  a  comprehensive  picture  of  the  displacement  mechanism  of 
the  oil  by  water  in  the  presence  of  the  free  gas.  The 
pressure  build-up  associated  with  the  water  flood  may 
completely  dissipate  the  free  gas,  reducing  it  to  a  two 
phase  system.  The  picture  may  also  be  complicated  by  the 
simultaneous  effects  of  partial  solubility  and  compressibility 
of  the  free  gas. 
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NO  GAS  PRESENT  GAS  PRESENT 

CZ3  SAND 
&Z3  WATER 
□  GAS 
WM  OIL 

FIG.  I.  CONCEPT  OF  OIL  AND  GAS  LOCATIONS  IN  PORE 
SPACE  AT  RESIDUAL  OIL  SATURATION. 

(FROM  REFERENCE  2.) 
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In  view  of  the  foregoing  factors,  the  investigation 
was  conducted  at  relativelyalow  (0.156  psi/cm.)  pressure 
differential  to  reduce  the  effects  of  solubility  and 
compressibility. 
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LITERATURE  REVIEW 


Laboratory  investigations  have  shown  that  there  is 
generally  an  increase  in  the  productions  of  oil  by  the  main¬ 
tenance  of  free  gas  saturations  during  the  flooding  opera¬ 
tions.  Investigators  have  presented  relationships  between 
initial  gas  saturation  and  oil  recovery,  as  well  as  between 
initial  gas  saturation  and  final  gas  saturation.  They,  1, 

2,  8,  10,  15,  etc.,  have  also  tried  to  calculate  the  recovery 

by  water  flood  with  the  use  of  relative  permeability  data. 

% 

Holmgrem  and  Morse  (2)  were  among  the  first  to  show 
the  relationship  between  the  initial  free  gas  saturation  and 
the  residual  oil  saturation.  They  also  presented  a  relation¬ 
ship  between  the  initial  free  gas  saturation  and  the  trapped 
gas  saturation  after  the  oil  floods.  In  their  work  the  free 
gas  saturation  was  created  by  flowing  gas  through  the  system. 
Once  the  desired  gas  saturation  was  obtained,  gas  saturated 
oil  was  injected  into  the  core  under  a  low  pressure  differential 
in  order  to  displace  excess  free  gas  thus  preparing  the  core 
for  water  flooding.  The  relationships  which  they  developed 
are  shown  in  Figures  2  and  3. 

Geffen  and  others  (20)  concluded  from  their  studies 
that  a  substantial  amount  of  gas  would  remain  in  the  pore 
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0  10  20  30  40  50 

TRAPPED  GAS  SATURATION  AFTER  OIL 

FLOOD, 

Percent  pore  volume. 


FIG.  2.  RELATIONSHIP  BETWEEN  MAXIMUM  GAS  SATURATION 
OBTAINED  BY  GAS  DRIVE  AND  TRAPPED  GAS  SATU¬ 
RATION  AFTER  OIL  FLOOD. 


(FROM  REFERENCE  2) 


7 


FIG.  3.  RELATION  BETWEEN  FREE  GAS  SATURATION  AND 
RESIDUAL  OIL  SATURATION  AT  IOO- 1  WATER  OIL 
RATIO 

CURVE  A  -  GAS  TRAPPED  BY  OIL  FLOOD  ONLY; 
CURVE  B  -  GAS  TRAPPED  BY  OIL  FLOOD  AND  EX¬ 
PANDED  BY  PRESSURE  REDUCTION, 

POINT  C-  GAS  ESTABLISHED  BY  INTERNAL  DRIVE. 


(FROM  REFERENCE  2.) 
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space  of  water- invaded  media  originally  containing  a  free 
gas  saturation.  It  was  estimated  that  the  range  of  trapped 
gas  saturation  would  be  15  to  50  percent  of  the  pore  space 
and  not  the  critical  gas  saturation  usually  taken  as  1  to 
11  percent  pore  space. 

The  purpose  of  a  study  conducted  by  Dyes (21)  was  to 
show  that  the  production  of  the  water  driven  reservoirs  at 
pressures  below  the  bubble  point  would  improve  the  oil  re¬ 
placement  efficiency  as  a  result  of  gas  evolution.  Figure 
4  gives  the  relationship  between  initial  and  final  gas 
saturations  and  Figure  5  presents  the  effect  of  gas  saturation 
on  the  recovery  of  oil  showing  maximum  recovery  for  a  definite 
range  of  pressure.  The  experiments  were  conducted  on  small 

cores  with  high  flooding  rates. 

Kennedy,  Kruger  and  others  (11)  investigated  the 
relationships  between  residual  gas  saturation  after  the  water 
floods  and  initial  gas  saturation  and  arrived  at  similar 
results.  They  also  studied  the  effects  of  saturation  history. 

Kyte  and  others  (6)  conducted  experiments  covering  a 
wide  range  of  core  materials  and  fluid  properties.  They 
investigated  the  mechanism  of  oil  displacement  by  water  in 
partially  gas- saturated  porous  media.  They  found  that  a 
high  viscosity  oil  gave  lower  residual  gas  saturation  while 
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TRAPPED  GAS  SATURATION  AFTER 
WATER  FLOOD, 

Percent  pore  volume. 


FIG-  4.  RELATIONSHIP  BETWEEN  INITIAL  AND  RESIDUAL 
GAS  SATURATION  FOR  TWO  SANDSTONE  CORE 
GROUPS. 


(FROM  REFERENCE  21.) 


OIL  RECOVERY  BY  SOLUTION  GAS  DRIVE  PLUS  WATER 
DISPLACEMENT,  Percent  initial  oil  in  place. 
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FIG.  5.  OIL  RECOVERY  BY  PRODUCTION  UNDER  SOLUTION 

GAS  DRIVE  FOLLOWED  BY  DISPLACEMENT  BY  WATER. 
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a  low  viscosity  oil  gave  a  higher  residual  gas  saturation. 

These  experiments  were  conducted  on  Torpedo  sandstone  and 
the  gas  used  was  helium.  It  was  also  concluded  in  this 
investigation  that  under  similar  initial  saturations  and 
flooding  conditions  the  residual  gas  saturation  after  a  water 
flood  would  be  less  for  an  oil-wet  system  than  for  a  water- 
wet  system.  Some  of  these  results  are  illustrated  in  Figures 
6  and  7.  From  the  investigations  they  showed  that  for  one 
given  system  there  is  a  unique  relationship  between  the  trapped 
gas  saturation  and  the  initial  mobile  gas  saturation.  They 
also  showed  that  for  the  given  systems  there  was  a  unique 
relationship  between  the  oil  saturation  attained,  at  any 
flooding  stage  and  the  trapped  gas  saturation  existing  in  the 
porous  medium  during  the  flood.  Furthermore^  it  was  found 
that  there  was  a  unique  relationship  between  the  oil  satura¬ 
tion  attainable  at  any  one  flooding  stage  and  the  initial 
mobile  gas  saturation. 

Dardaganian  (18)  conducted  experimental  work  and 
obtained  similar  results  regarding  the  relationships  between 
initial  and  final  gas  saturation.  He  also  investigated  the 
effects  of  viscosity  ratios  and  the  parameter  k/o.  The 
results  are  illustrated  by  Figures  8,  9,  and  10. 
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FIG.  6.  DISPLACEMENT  OF  MOBILE  GAS  BY  OIL  AND  WATER 
INJECTION. 


(FROM  REFERENCE  6.) 


OIL  SATURATION,  Percent  pore  volume. 
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(FROM  REFERENCE  6.) 


VISCOSITY  OF  THE  DISPLACING  PHASE, 
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INITIAL  GAS  SATURATION, 
Percent  pore  volume- 


FIG-  8-  EFFECT  OF  VISCOSITY  RATIO  (p0/pg)  0N 
RESIDUAL  GAS  SATURATION- 


(FROM 


REFERENCE  18) 


RESIDUAL  GAS  SATURATION, 
Percent  pore  space. 
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INITIAL  GAS  SATURATION, 

Percent  pore  volume- 

FIG-  9-  RELATIONSHIP  BETWEEN  RESIDUAL  AND  INITIAL  GAS 
SATURATION  FOR  SYSTEMS  OF  VARYING  ROCK 

PARAMETER  (K/0f).* 
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Bass  (7)  concluded  from  his  experimental  work  that  a 
peak  recovery  is  obtained  by  water  flooding  at  an  initial 
gas  saturation  of  5  to  7  percent  and  the  pressure  at  which 
oil  recovery  is  obtained  is  insensitive  to  initial  water 
saturations  within  the  range  investigated. 

Among  the  investigators  who  presented  methods  for 
the  calculation  of  recoveries  by  water  flooding  in  the 
presence  of  free  gas  are  Holmgren  and  Morse  (2),  Pirson  (5), 
Stulzman  and  Thodos  (14),  Donglas  Jr.  (15),  Welge  (8), 

Rapoport  and  Lea  (9)  and  Dardagamian  (10). 

Holmgren  and  Morse  (2)  calculated  water  flood 
recovery  by  using  Buckley  and  Leveretts  (1)  method  as 
simplified  by  Pirson(5) .  For  the  calculation  of  recovery 
at  breakthrough  for  a  trial  run  they  used  relative  permeability 
to  oil  and  water  for  a  trapped  gas  saturation  of  25  percent. 
Their  trial  calculation  was  in  close  agreement  with  the 
experimental  result. 

Donglas  Jr.  and  others  (15)  presented  a  numerical 
method  to  determine  the  behavior  of  a  linear  flood  with  the 
inclusion  of  the  effects  of  capillary  pressure.  Leverett 
(19)  originally  derived  a  second-order,  non  linear  duferen- 
tial  equation  including  the  capillary  pressure  term  for  the 
displacement  of  oil  from  a  linear  porous  medium  by  an 
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immiscible  fluid.  This  equation  was  approximated  by  an 
implicit  form  of  the  difference  equation  for  solution  by 
the  electronic  digital  computor.  From  their  experimental 
work,  using  an  oil-wet  system  they  showed  that  the  efficiency 
of  the  flood  and  the  sharpness  of  the  flood  front  initially 
increased  with  an  increasing  rate  and  tnen  levellea  off 
with  no  further  change.  However,  for  a  water -wet  system, 
the  breakthrough  recovery  passed  through  a  minimum  and  then 

increased  with  decreasing  rate  of  flood. 

Welge  (8),  starting  with  the  fundamental  relationships 
given  by  Buckley  and  Leverett  (1),  derived  a  useful  analyti¬ 
cal  method  for  computing  the  breakthrough  recovery.  The 
advantage  of  the  method  is  in  the  use  of  relative  permeability 
data  in  a  limited  range.  The  tangential  construction  explained 
in  Figure  14  by  Welge  presents  a  simplified  method  for  the 
calculation  of  time  of  breakthrough  and  breakthrough  recovery. 

Dardaganian  (10)  presented  a  simple  mathematical 
proof  for  the  basic  tangential  construction  presented  by 
Welge  (8)  and,  in  addition,  modified  it  for  various  initial 
saturation  conditions.  Buckley  and  Leverett' s  fractional 
flow- frontal  advance  theory  breaks  down  when  applied  to  a 
system  with  three  mobile  phases.  Dardaganian  broke  the  three 
phase  system  down  to  two  separate  systems  in  which  only  two 
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mobile  phases  exist.  One  system  is  that  in  which  the  oil 
bank  displaced  the  gas  thereby  reducing  it  to  an  immobile 
saturation  in  the  presence  of  immobile  water  phase.  In 
the  second  system  the  water  flood  front  displaced  the  oil 
from  the  porous  rock  in  the  presence  of  an  immobile  gas 
phase.  The  residual  oil  saturation  by  water  aisplacement 
is  a  function  of  the  residual  gas  within  the  oil  bank.  The 
effect  of  the  residual  gas  saturation  on  the  resiaual  oil 
was  computed  by  the  fractional  flow  -  frontal  advance  theory 
with  appropriate  relative  permeability  (ko/kw)  data  for  various 


free  gas  saturations. 
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EXPERIMENTAL  EQUIPMENT 

A  schematic  diagram  of  the  experimental  apparatus 
is  presented  in  Figure  11.  A  cylindrical  core  of  Berea 
sandstone  32  centimeters  long  and  5  centimeters  in  diameter, 
mounted  in  lucite  and  equipped  with  eight  pressure  taps, 
was  the  test  sample.  A  duplex  Milton  Roy  pump  was  used  to 
pump  the  live  oil  from  a  seven  litre  stainless  steel  cylinder 
into  the  system.  Glass  wool  filters  were  fitted  in  the  lines 
at  the  inflow  end  of  the  core.  A  pressure  regulator  at  the 
methane  bottle  kept  a  constant  30  psig  gas  pressure  on  the 
live  oil  which  consisted  of  -72°F  pour  point  diesel  saturated 
with  99.8  percent  pure  methane.  One  120  cc.  graduated  lucite 
tube  was  used  to  hold  the  flooding  brine.  A  Grove  back- 
pressure  regulator  fitted  at  the  effluent  end  of  the  core 
was  used  to  control  all  the  fluid  transfers  and  pressures 
in  the  core.  A  lucite  separator  fitted  with  thirteen  15  ml. 
graduated  test  tubes  and  two  100  ml.  gas  measuring  burettes 
which  were  fitted  with  three-way  valves  were  used  to  measure 
the  effluents  at  atmospheric  pressure. 
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A  temperature  regulated  cabinet  held  all  the  equipment 
except  the  Milton  Roy  pump,  the  methane  bottle  and  the 
effluent  measuring  devices. 
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EXPERIMENTAL  PROCEDURE 

Ten  thousand  parts  per  million  sodium  chloride  solu¬ 
tion  in  water  was  used  for  connate  water  as  well  as  for  the 
flooding  fluid.  Diesel  oil  with  a  -72°F  pour  point  was 

saturated  with  pure  methane  at  30  psig. 

To  prepare  the  core  for  the  reference  run,  i.e.  for 
zero  percent  initial  gas  saturation,  the  core  was  evacuated 
and  completely  saturated  with  brine  under  pressure.  The 
brine  was  then  displaced  from  the  core  with  live  oil  by 
pumping  it  through  the  core  at  a  pressure  well  above  the 
bubble  point.  When  the  production  ratio  of  oil  to  water  was 
above  1000,  the  pump  was  stopped  and  the  core  was  ready  for 
the  reference  run.  The  connate  water  saturation  of  the  coie 

thus  obtained  was  49.5  percent. 

All  the  water  floods  were  conducted  at  five  pounds 
pressure  differential  across  the  core.  An  air  pressure  source 
was  used  to  inject  the  brine.  The  back  pressure  was  always 
set  equal  to  the  average  pressure  in  the  core.  The  volumes 
of  the  brine  injected  and  the  effluents  were  recorded  during 
the  flood  and  also  during  the  preparatory  steps.  Furthermore, 
at  the  inlet  and  outlet  and  in  the  sections  of  the 


pressures 


1 


■ 


24 


core  were  recorded  as  a  function  of  time  during  the  water 
floods.  Flooding  was  stopped  at  a  water-oil  ratio  of  100. 

In  preparation  for  runs  with  initial  gas  saturation 
in  the  system,  the  core  at  the  end  of  a  water  flood  was 
flooded  with  live  oil  at  a  pressure  well  above  the  bubble 
point.  The  oil  flood  was  stopped  when  the  predetermined 
connate  water  saturation,  known  from  a  material  balance,  was 
reached. 

Before  each  run  the  initial  pressure  in  the  core 
saturated  with  live  oil  was  50  psig.  The  core  was  depleted 
at  an  arbitrary  rate  of  0.25  ml. /min.  until  sufficient  liquid 
had  been  produced  to  create  the  required  gas  saturation. 

The  core  was  then  left  for  18  hours  for  stabilization  and  to 
minimize  the  effects  of  super saturation. 

To  create  gas  saturation  above  that  possible  by  de¬ 
pletion  drive,  gas  drive  was  conducted  at  the  end  of  depletion 
drive.  Keeping  the  material  balance  on  liquids  the  predeter¬ 
mined  gas  saturation  was  thus  obtained. 

^11  the  floods  were  conducted  at  100°F  and  the  effluents 

were  collected  at  74°F  and  13.7  psia. 
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EXPERIMENTAL  RESULTS  AND  DISCUSSION 

The  data  obtained  in  this  work  are  shown  in  graphical 
form  and  are  tabulated  in  the  Appendix.  The  results  are 

presented  in  Tables  1  and  2. 

A  theoretical  calculation  of  the  reference  run  (water 
flood  at  zero  gas  saturation)  was  made  to  check  the  experimen¬ 
tal  results. 

Figure  15  presents  a  tangential  construction  required 
for  the  calculation.  The  viscosity  ratio  of  oil  to  water  in 
the  reference  run  was  2.65.  The  tangential  construction  is 
fully  explained  in  the  Appendix.  This  construction  permits 
determination  of  the  average  water  saturation  at  breakthrough 
as  well  as  the  water  saturation  at  the  flood  front.  Relative 
permeability  data  used  were  found  by  pumping  oil  above  its 
bubble  point  pressure  through  the  core  saturated  with  oil 
and  brine.  The  breakthrough  recovery  obtained  graphically 
was  19.6  percent  of  the  pore  volume  as  compared  to  21.7  per¬ 
cent  obtained  experimentally.  The  difference  in  the  values 
of  oil  recovery  obtained  graphically  and  experimentally  is 
probably  a  result  of  unsatisfactory  relative  permeability 
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Run  No.  10  includes  depletion  drive  plus  gas  drive. 

Rate  of  Production  =  0.25  ml. /min. 

Time  of  Stabilisation  =  18  hours. 

Correction  applied  for  the  production  above  bubble  point  pressure 
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data  and  the  rate  sensitivity  of  the  flood.  The  latter 
will  be  discussed  later. 

Figure  16  shew  s  the  relationship  between  gas  satura¬ 
tion  at  the  start  of  the  flood  and  the  residual  gas  saturation 
at  the  end  of  the  flood.  Curve  A  shows  the  above  mentioned 
relationship  with  the  residual  gas  saturation  as  it  existed 
in  the  rock  at  the  average  pressure  at  the  end  of  the  water 
flood  and  curve  B  is  the  residual  gas  saturation  at  the  average 
pressure  of  the  core  at  the  start  of  the  water  flood.  From 
curve  B  we  see  that  the  residual  gas  saturation  is  equal  to 
the  initial  gas  saturation  from  zero  to  almost  10  percent  gas 
saturation.  Above  that  there  is  a  proportional  increase  of 
residual  gas  saturation  with  the  increase  of  initial  gas 
saturation.  For  an  initial  gas  saturation  of  33.4  percent 
pore  volume  the  residual  gas  saturation  is  20.7  percent  pore 
volume.  This  phenomenon  is  in  close  agreement  with  the 
results  obtained  by  Holmgren  and  Morse  (2)  and  Kruger  and 
others  (11)  but  does  not  completely  agree  with  Dyes (21)  in 
the  lower  range  of  gas  saturation.  In  the  very  low  range  of 
gas  saturation  Dyes (21)  does  not  find  the  residual  gas  satura¬ 
tion  equal  to  the  initial  gas  saturation  although  it  is  very 
close.  This  difference  probably  is  a  result  of  the  charac¬ 
teristics  of  the  rock.  In  some  cases  part  of  the  gas  createc 
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FIG.  16.  RELATIONSHIP  BETWEEN  RESIDUAL  AND  INITIAL  GAS 

SATURATION. 
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GAS  SATURATION  AT  START  OF  WATER  FLOOD,  Percent  pore  volume. 
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in  the  rock,  however  little  it  may  be,  is  displaced  with 
the  movement  of  the  oil  while  in  other  cases  it  is  not. 

Figure  17  presents  plots  of  oil  saturations  in  the 
core  during  the  flood  versus  the  cumulative  brine  injected. 

It  can  be  noticed  that  curves  0,  1,  2  and  3  are  almost 
parallel  to  each  other.  The  trend  starts  changing  with 
curves  4  and  5.  Curves  6  to  10  have  definite  horizontal 
portions  which  increase  with  the  increase  of  the  initial 
gas  saturation.  For  curves  0  to  4  the  production  of  oil 
starts  almost  instantaneously  with  the  invasion  of  water, 
while  for  higher  gas  saturations  there  is  a  time  lag  between 
the  initiation  of  the  flood  and  the  production  of  oil.  This 
time  lag  increases  with  increasing  initial  gas  saturation  as 
shown  by  the  flat  portions  of  the  curve.  It  might  be  con¬ 
cluded  from  these  results  that  during  water  flooding  two 
types  of  displacement  mechanism  exist  which  are  governed  by 
the  saturation  of  the  gas  phase  present  at  the  start  of  the 
flood.  For  lower  gas  saturations  the  basic  mechanism  of 
displacement  of  oil  by  water  is  the  same  as  for  one  without 
any  gas  phase.  The  gas  present  in  the  system  is  immediately 
trapped  by  the  movement  of  the  oil  and  the  system  reduces  to 
water  displacing  oil  in  the  presence  of  immobile  water  and 
gas  phases  instead  of  immobile  water  alone.  The  only  effect 


FIG-  17.  RELATIONSHIP  BETWEEN  BRINE  INJECTED  AND 

OIL  SATURATION. 
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of  trapped  gas  would  be  to  decrease  the  relative  permeability 
of  oil.  For  higher  gas  saturations  the  mechanism  changes  to 
a  system  which  can  be  separated  into  two  distinct  processes: 
the  oil  bank  displacing  the  gas  phase  leaving  behind  a 
trapped  gas  in  the  presence  of  immobile  water  followed  by 
water  displacing  oil  in  the  presence  of  trapped  immobile  gas. 
This  change  in  the  mechanism  for  the  system  takes  place  at 
8  to  10  percent  pore  volume  gas  saturation  for  the  system 
studied. 

It  is  usually  visualized  that  during  a  water  flood 
most  of  the  gas  present  in  the  system  goes  back  into  solution 
behind  the  oil  bank  and  in  front  of  the  water  bank.  From 
the  material  balance  calculations  on  the  liquid  phases  it 
is  seen  that,  at  such  low  pressure  differentials  as  encountered 
in  this  experimental  work,  the  amount  of  gas  going  back  into 
solution  is  insignificant.  The  reason  probably  is  the  lack 
of  time.  Also,  the  volume  change  of  the  oil  due  to  the 
solubility  of  the  gas  is  insignificant  as  can  be  observed 

from  Table  7  in  the  Appendix. 

Figure  18  presents  a  relationship  between  brine  in¬ 
jected  and  the  oil  production.  Curves  0  to  6  are  close  to¬ 
gether  and  all  converge  towards  the  origin  while  curves  7  to 
10  show  horizontal  parts  which  increase  with  the  increase 
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of  gas  saturations  at  the  start  of  the  water  flood.  These 
results  (shown  in  Figure  18)  substantiate  the  discussion  of 
the  two  different  mechanisms  of  the  water  flood. 

Figure  19  presents  the  relationship  between  brine 
injected  and  gas  saturation.  Curves  1  to  4  are  almost 
horizontal  and  parallel  lines  showing  no  effect  of  change 
of  gas  saturation  on  the  mechanism  of  displacement.  The 
trend  of  the  curves  from  5  to  10  changes  with  initial  gas 
saturation . 

Figures  21a  to  21 j  present  the  pressure  profiles  of 
the  core  taken  during  the  water  floods.  A  specimen  pressure 
profile  for  water  flood  number  0  has  been  presented  in  this 
section  in  Figure  12.  The  arrows  represent  the  possible 
location  of  the  water  flood  front.  The  advance  of  the  increased 
pressure  drop  zone  is  indicative  of  the  movement  of  the  flood 
front.  This  pressure  drop  is  caused  by  the  sudden  change  of 
capillary  forces  and  relative  permeabilities  at  the  flood 
front.  Since  the  pressure  differential  of  five  pounds  across 
the  core  was  quite  small,  it  was  sometimes  difficult  to  lo¬ 
cate  the  flood  front  from  the  pressure  profiles,  especially 

at  higher  gas  saturations. 

The  most  important  relationship  of  initial  gas 
saturation  and  the  total  recovery  (depletion  drive  plus 
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FIG.  12.  PRESSURE  PROFILES  DURING  WATER  FLOOD  No.  0. 
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water  flood)  is  presented  in  Figure  20.  A  recovery  of 
42.75  percent  of  initial  oil  in  place  is  obtained  in  the 
reference  water  flood  (in  the  absence  of  any  gas).  The 
recovery  of  the  oil  increases  with  the  increase  of  gas 
saturation  before  the  water  flood.  At  the  initial  gas 
saturation  of  33.4  percent  pore  volume  the  total  recovery  is 
64.65  percent  of  initial  oil  in  place. 

This  result  does  not  agree  with  Bass  (7),  Dye  (21), 
Holmgren  or  Dardaganian  (18)  though  similar  results  were 
obtained  by  Saxon  and  others  (3).  Bass  (7)  deduces  from 
his  studies  that  under  the  conditions  of  his  investigation 
a  peak,  oil  recovery  is  obtained  by  water  flooding  at  initial 
gas  saturation  of  5  to  7  percent  pore  volume. 

The  fractional  flow  equation  developed  by  M.C. 
Leverett  (19)  is 


1  + 

It 


n  (6PC/6U)  -  g  A P  Sin© 
pQ  


1  + 


\h £ 

P'0 


(A) 


where  fw  -  fraction  of  the  displacing  fluid  in  the  flow 

stream 

K  -  effective  permeability 

p,  -  viscosity 

P  -  density 

q  -  total  fluid  rate  per  unit  cross-sectional  area 


FIG.  20.  RELATIONSHIP  BETWEEN  GAS  SATURATION  BEFORE 
WATER  FLOOD  AND  TOTAL  RECOVERY 


38 


o 

cn 


30D|d  Ul  |I0  |Dj|IU|  |U90J8d 


lAd3A0D3d  3AIHQ  d31VM  Sffld  Ad3A003d  3AldQ  NOI13"ld3Q 


GAS  SATURATION  AT  START  OF  WATER  FLOOD, 
Percent  pore  volume. 


39 


g 


u 


P  =P  -P 
rc  o  rw 

0 


-  gravity 

-  length  in  the  direction  of  flow 

-  pressure  drop  in  the  two  fluids 

-  dip  angle  from  the  horizontal 


Subscript  o  -  oil 


w  -  water 


From  this  it  can  be  seen  that  fractional  flow  of  the 
displacing  fluid  is  a  function  of  rate  and  capillary  forces 
as  well  as  other  factors.  In  addition,  consideration  of 
the  dimensionless  equations  (B  and  C)  developed  by  Rapoport 
and  Leas  (9)  indicates  that  the  behavior  of  the  linear  flood 
is  a  function  of  the  length  of  the  flooded  system  and  the 


rate  of  injection. 

ds  +  dF  5s  k 

"St  ds  ’  5x 


LVp 


w 


5 

3x 


K0  .  F  . 


dP( 

ds 


5s 

3x 


=  o  (B) 


K, 


1  + 


LVpw 


^w 


dP( 

ds 


5s 

3x 


=  0  at  X  =  0,  for  any  T  (C) 


where  k  -  specific  permeability 

_  relative  permeability  to  water  and  oil  respectively 

p0  -  viscosities  of  water  and  oil  respectively 
s  -  water  saturation 
f  -  porosity 
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V  -  total  flow  rate  per  unit  cross-section  area 


c  - 
F  - 
X  - 
L  - 


M-o  /  f^w 


K, 


-1 


1  +/ — £JX_\ 

(  C  V 
x/L,  x  -  distance 

total  length  of  flooded  system 


T  -  tV/Lf 
t  -  time 

LVpw  -  scaling  coefficient 
Pc  -  capillary  pressure 

Accordingly  all  floods  corresponding  to  the  same  scaling 
factor,  which  is  the  product  of  length  of  the  system,  velocity 
of  the  flood  and  viscosity  of  the  displacing  fluid,  must 
behave  in  a  similar  manner  and  yield  equal  recoveries  for 
equal  cumulative  injection.  The  experiments  were  conducted 
on  an  oil-wet  alundum  core.  Donglas  Jr „  05)  also  investigated 
the  effect  of  rate  of  flood  on  the  recovery  for  homogeneous 
systems.  According  to  him,  in  a  water  wet  system,  the  break¬ 
through  recovery  increases  with  decreasing  rate  of  water  flood 
for  low  flooding  rates.  Such  a  behavior  is  a  result  of  the 
spreading  of  the  invading  water  front  by  capillary  forces  and 
the  end  effect.  At  extremely  low  rates  the  latter  effect 
would  predominate  and  prevent  the  production  of  water  until 
all  the  recoverable  oil  has  been  produced. 
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Since  the  creation  of  gas  in  the  system  seems  merely 
a  replacement  of  oil  by  gas  in  the  oil  phase  it  is  logical 
to  conclude  that  breakthrough  recovery  will  increase  with 
an  increase  of  the  gas  phase  in  the  system  for  low  rates 
where  capillary  forces  and  end  effects  predominate.  Such 
is  the  case  for  the  system  investigated. 


CONCLUSIONS 


Under  the  conditions  of  this  work  the  following 
conclusions  may  be  drawn: 

1.  An  increase  in  total  recovery  is  obtained  by  water 
flooding  in  the  presence  of  free  gas  in  the  system  rather 
than  by  flooding  in  the  absence  of  the  gas  phase. 

2.  Up  to  33.4  percent  gas  saturation  the  total  recovery  of 
oil  by  water  flooding  increases  with  an  increase  of  initial 
gas  saturation.  A  definite  increase  of  approximately  22 
percent  of  initial  oil  in  place  is  obtained  by  increasing 
initial  gas  saturation  from  zero  to  33.4  percent. 

3.  An  initial  gas  saturation  up  to  approximately  10  percent 
does  not  basically  change  the  mechanism  of  displacement. 

4.  In  comparing  with  other  results  it  may  be  concluded  that 
for  each  separate  system  there  is  a  unique  relationship 
between  initial  and  trapped  gas  saturation. 

5.  A  gas  drive  may  be  conducted  before  the  water  flood  so 
as  to  obtain  higher  gas  saturations. 
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APPENDIX 


The  fundamental  development  of  the  Buckley-Leverett 

theory  by  Dardaganian  (10)  is  as  follows. 

Neglecting  effects  due  to  gravity  and  capillary 

pressure  differences  the  factional  flow  equation  originally 

developed  by  M.C.  Leverett  (19)  may  be  expressed  as 

* 

1 


fd  = 


ko  iid 
1  +  |io  kd 


(1) 


The  frontal  advance  relationship  developed  by  Buckley 
and  Leverett  (1)  is  expressed  in  equation  2 


Ax  =  Qt  Afd, 
AT  Kb 


(2) 


The  basic  mechanism  of  fluid  displacement  has  been 
discussed  at  many  placed  in  the  literature.  Qualitatively 
this  mechanism  may  be  described  by  the  saturation  distribution 
represented  in  Figure  13.  The  abscissa  is  the  pore  volume 
and  the  ordinate  is  the  displacing  phase  saturation  in  percent 

pore  space. 


*  Nomenclature  is  presented  at  the  end  of  this  section 
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For  the  purpose  of  discussion  the  system  presented 
was  considered  to  contain  two  phases,  oil  and  the  displacing 
phase.  Prior  to  injection,  the  saturation  in  the  displacing 
phase  was  considered  to  be  Sdi;  therefore,  the  initial  oil 
saturation  was  1  -  Sdi.  The  oil  saturation  at  which  the 
permeability  to  oil  vanished  in  a  relative  permeability  test 

was  considered  to  be  Sor. 

The  displacing  phase  accumulates  in  the  system  re¬ 
placing  the  oil.  As  the  displacing  fluid  passes  into  the 
system,  a  front  or  stabilized  zone  of  fixed  saturation  dis¬ 
tribution  is  developed  within  which  each  plane  of  saturation 
moves  with  a  velocity  dependent  on  the  value  of  saturation. 

At  the  inflow  face  of  the  system  the  oil  is  reduced  to  a 
value  of  Sor  since  many  pore  volumes  of  displacing  fluid  pass 

the  inflow  face  in  a  short  time. 

If  both  phases  are  mobile  at  initial  conditions  and 
considering  the  system  to  be  in  a  steady  state,  then  at  the. 
out  flow  face  Qt  =  Qtfoi  +  Qt  fdi  until  breakthrough.  That 
is,  the  composition  of  the  flow  stream  from  the  system  is 
determined  by  the  saturation  prevailing  at  the  outflow  end 

of  the  system. 

At  a  given  time  T,  the  volume  of  the  fluid  injected 
into  the  system  is  given  by  QtT;  therefore,  the  volume  of 
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displacing  fluid  accumulating  in  the  system  is  equal  to 
QtTfoi-  Designating  the  average  saturation  in  the  dis¬ 
placing  phase  behind  the  front  as  Sda,  then  by  vulume 
balance,  if  the  front  has  advanced  to  position  x  at  the 


end  of  T  units  of  time, 

QtTfoi  =  x  AO  (Sda  -  Sdi) 

From  Figure  13  it  may  be  noted  that 


(3) 


QtTfoi  =  x  M  (Sdf  “  Sd0  + 

J  Suf 


1-Sor 


x  A<t>  dSd  (4) 


Also  from  Equation  2 


x  '  W  Sd 


and 


5  -  *£  <S> 


sdf 


(3) 


(6) 


Substituting  equation  (5)  into  equation  (4)  and  integrating 
QtT  foi  =  x  A<J>  (Sd£  -  Sdi)  +  QtT  (1  "  fdf) 

di 


Since  fo£  =  1  -  fdi 


QtT  "  fd0  =  *  A0  (Sdf 
From  equation  (6) 

✓ 


-  sdi)  +  QtT  a  -  fdf)  (8> 


Q  T  =  x  A$ 


dfd 

(dSi} 


(9) 


»df 
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Therefore 
, _ 1 

x  A<i> 


df  d 
( — -) 
W 


(i-fdi)  =  x  A4>  (Sdf-Sdi)  +  x  A4> 


Sdf 


dfd 


d  S 


0--fdf) 


df 


or 


(f!i) 

dSd  Sdf 


(l-fdi)  -  (l-fdf) 


and 

(f£d)  =  f4f  ' 

dsd  Sdf  -  sd^ 


(10) 


Also  combining  equations  (3)  and  (6) 


) 


Sdf 


or 


1  ~  fdi 
sda  “  Sdi 


£df  ~  £di 

sdf  "  sdi 


(ID 


(12) 


The  slope  of  the  fractional  flow  curve  at  the 
saturation  of  the  front  must  thus  satisfy  equation  10  and 
equation  11  which  state  that  the  tangent  line  must  pass 

through  the  points  (Sdi,  fdi^  (sdf  “  fdf)  and  (Sda-»  * 

The  saturation  at  the  front  and  the  average  saturation 

behind  the  front  may  thus  be  defined  by  drawing  a  fractional 
flow  curve  and  constructing  a  tangent  line  to  the  curve  with 
the  initial  point  having  the  coordinates  (Sdi,  f di) .  The 
saturation  corresponding  to  that  at  the  front  is  given  by 
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the  point  of  tangency  and  that  corresponding  to  the  average 
saturation  behind  the  front  by  the  intercept  with  the  ordinate, 
fd  =  1.  (See  Figure  14) 

It  may  further  be  noted  that  only  one  such  tangent 
line  can  be  constructed  from  a  given  set  of  initial  coordinates. 
Also  it  is  apparent  that  if  the  initial  saturation  in  the 
displacing  phase  exceeds  that  corresponding  to  the  maximum 
slope  of  the  fractional  flow  curve  no  tangent  line  may  be 
drawn  satisfying  the  stated  conditions.  This  is  interpreted 
to  mean  that  a  front  will  not  be  developed  in  the  system. 

The  oil  recovery  at  breakthrough  expressed  in  percent 
of  pore  space  is  given  by  (Sda  -  Sdi)  or,  m  percent  of  oil 

originally  in  place  by  S^a  -  _  Sda  -  Sd^  .  a  similar 

Soi  1  “  Sdi~ 

approach  for  the  problem  with  three  phases  has  been  discussed 
in  the  previous  sections. 

Figure  15  presents  a  similar  tangential  construction 
drawn  for  reference  water  flood. 
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DISPLACING  PHASE  SATURATION,  Sd 


TANGENTIAL  CONSTRUCTION  WITH  INITIAL  DISPLACING 
PHASE  SATURATION  IN  TWO  PHASE  SYSTEM. 


(FROM  REFERENCE  IO-) 
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FIG.  15-  FRACTIONAL  FLOW  CURVE  FOR  WATER  FLOOD  No.  0 
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NOMENCLATURE 

A  -  cross-sectional  area  normal  to  the  direction  of  flow 
4>  -  porosity 

f  ^  _  fraction  of  the  displacing  fluid  in  the  flowing  stream 
d 

k  -  relative  permeability 

Qt  -  total  flow  rate 
T  -  time 

x  -  distance  along  the  path  of  flow 

x  -  distance  along  the  path  of  flow  to  the  flood  front 

p  -  viscosity 
S  -  fluid  saturation 

Subscripts 
a  -  average 
d  -  displacing  phase 
f  -  front 
g  -  gas 
i  -  initial 
o  -  oil 
r  -  residual 
t  -  total 
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FIG.  21b.  WATER  FLOOD  No.  2. 
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FIG.  21c.  WATER  FLOOD  No.  3. 
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FIG.  2ld  WATER  FLOOD  No.  4. 
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FIG.  2le.  WATER  FLOOD  No.  5. 
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FIG.  2lf.  WATER  FLOOD  No.  6. 
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FIG.  2lg.  WATER  FLOOD  No.  7. 
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FIG.  2lh.  WATER  FLOOD  No.  8. 
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FIG.  2li.  WATER  FLOOD  No.  9. 
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FIG.  2lj.  WATER  FLOOD  No.  10- 
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FIG.  23.  RELATIONSHIP  BETWEEN  VISCOSITY  AND  PRESSURE. 
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TABLE  3 


GENERAL  PROPERTIES  OF 

THE 

CORE  * 

Length  of  the  core 

= 

32.0 

cms . 

Diameter  of  the  core 

— 

5.0 

cms . 

Weight  of  the  dry  core 

— 

5528.5 

gms. 

Pore  volume 

— 

143.5 

c„c. 

Porosity 

— 

22.8 

percent 

Liquid  permeability 

— 

156 

md. 

*  Berea  Sandstone 
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End  of  the  run 

Breakthrough  time  =  156  minutes 
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End  of  the  run 

Breakthrough  time  =  138  minutes 
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Breakthrough  time  =  115  minutes 
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Breakthrough  time  =  126  minutes 
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End  of  the  run 

Breakthrough  time  =  128  minutes 
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End  of  the  run 

Breakthrough  time  =  135  minutes 
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End  of  the  run 

Breakthrough  time  =  190  minutes 
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TABLE  6 

FRACTIONAL  FLOW  CALCULATIONS  FOR  THE  REFERENCE  RUN 


Number  1 

2 

3 

4 

5 

6 

7 

8 

Brine  Saturation, 

ml.  98.7 

95.7 

92.7 

89.7 

87.7 

86.7 

85.7 

84.1 

Brine  Saturation, 

7o  Pore  Volume  68.8 

66.7 

64.6 

62.5 

61.1 

60.4 

59.7 

58.6 

Rate  of  Production 
of  Brine,  ml . /sec . 0 . 025 

0.025 

0.0175 

0.0151 

0.0133 

0.010 

- 

0.0042 

Rate  of  Production 
of  Oil,  ml. /sec.  0.0 

0.0 

0.0075 

0.0099 

0.0117 

0.015 

- 

0.0208 

£  =  1/1+ko/ls,  . 

* 

- 

0.701 

0.655 

0.534 

0.400 

- 

0.167 

*  Viscosity  ratio  of  oil  to  water  for  the  run  2.65 
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TABLE  7 

PHYSICAL  PROPERTIES  OF  THE  FLUIDS 


Pressure,  Solubility*  of  Formation  Viscosity  of  Viscosity  of 
Psis  Methane  in  the  Volume  Diesel  plus  Methane  at 

Diesel,  cc./cc.  Factor  Methane  System  100  F,  cp. 


0.0 

0.0 

1.00 

9.4 

0.402 

1.0021 

14.7 

0.678 

1.0026 

21.3 

1.017 

1.0031 

27.9 

1.340 

1.0052 

41.9 

2.223 

1.0071 

Viscosity 

of  the  brine 

at  100°F 

Density  of 

the  brine  at 

100°F. 

Density  of  the  diesel  at  100°F. 


at  100  F,  cp. 


1.99 

0.0155 

1.925 

0.0156 

1.905 

0.01565 

1.875 

0.01575 

1.860 

0.01586 

1.830 

0.0161 

=  0.6947  cp. 

=  1.0  gm./ml. 

=  0.7865  gm./ml. 


*  Solubility,  cc.  at  13.7  psia  and  74  F./cc.  at  13.7  poia  and 
100°F. 


